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The conventional Clausius-Mossotti polarization equation of state is known to be unstable for 
polar liquids having molecules with high polarizability. Room temperature water is an important 
example. The instability in the polarization equation of state is of the typical loop form requiring an 
"equal area" construction for studying the stable ordered phase. The ordered phase of a Clausius- 
Mossotti polar liquid then consists of domains each having a net polarization. The polarization may 
vary in direction from domain to domain. The ordered phases are quite similar to those previously 
discussed on the basis of Dicke superradiance. 
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I. INTRODUCTION 

The Clausius-Mossotti[ll polarization equation of 
state is quite often discussed in text books on electro- 
magnetic theory 0, 0, 0. The discussion is thought to 
give students a good physical idea about how microscopic 
models of dipole moments yield non-trivial macroscopic 
material dielectric constants e. It is only sometimes men- 
tioned that the Clausius-Mossotti model exhibits an in- 
stability hinting at a phase transition into an ordered 
polarized state. The unstable precursor to the phase 
transition is in fact self evident. The instability is par- 
ticularly strong in polar liquids such as water. In fact, 
the Clausius-Mossotti model in its classic text book form 
indicates that water at room temperature and pressure 
is in an ordered polarization state. 

The ordered state of polar liquids is described in the 
literature as exhibiting superradiance]^, 0, E El • I n the 
superradiant state, the polarization P has a mean value 
roughly constant in magnitude within an ordered do- 
main. However, P may vary in direction as one goes 
from one ordered fluid domain to another ordered fluid 
domain. The theoretical derivation of superradiant do- 
mains normally starts from condensed matter quantum 
electrodynamics. Yet many insights into the nature of 
superradiant thermal equilibrium can be obtained from 
a classical electrodynamic viewpoint. Our purpose is to 
derive the polar liquid phase transition implicit in the 
Clausius-Mossotti model. Our arguments employ (apart 
from classical electrodynamics) the elementary statisti- 
cal thermodynamics previously worked out by Debye and 
Langevin|10l Hl| . 

In Sec. II, we briefly discuss the derivation of the 
Clausius-Mossotti local electric field (in Gaussian units) 



F = E+ (4tt/3)P. 



(1) 



In Eq.Q), E is the total electric field and P is the dipole 
moment per unit volume. In terms of the molecular po- 
larizability a of a single molecule and the liquid volume 
v per molecule, wc shall derive the mean field Clausius- 
Mossotti dielectric constant e. In Sec. Ill, we discuss 
the region of thermodynamic stability of the Clausius- 



Mossotti model. The known instability is often thought 
to be an embarrassment for the mean field theory. How- 
ever, we argue that the Clausius-Mossotti instability is a 
signature of a physical phase transition into an ordered 
phase. The ordered phase has been previously discovered 
from the viewpoint of superradiance. The thermody- 
namic theories of both the ordered and disordered phases 
are discussed from the Clausius-Mossotti viewpoint in 
Sec. IV. The nature of a more complete quantum theory 
is explored in the concluding Sec. V. 



II. CLAUSIUS-MOSSOTTI LOCAL FIELDS 

Consider the Coulomb energy contained in a charge 
density p, i.e. 



i r rp(^W) 

2./ ./ r-r' 



(2) 



If the charge density is derived from a macroscopic po- 
larization 

p(r) = -di«P(r), (3) 
then integrating Eq.J2J) twice by parts yields 

U = ~ j f P(r) • G(r - r') • P(r') d 3 rdV. (4) 



The dyadic 



G(r — r') = grad grad' 



r — r 



(5) 



obeys the trace condition 



TYG(r - r') = -V 2 r = Air8{r - r'). (6) 

r — r' 



Eqs.JSJ) and ijfjj imply that 



47T 



G(R) = — A(R) + — 1<5(R), 



A(R) 



3RR- R 2 1 



(7) 
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where the (47rl/3)<5(R) term in Eq.JJJ is present to en- 
force the trace condition in Eq.©. From Eqs.Q and 
0, it follows that the Coulomb energy 



U U dipole 4~ Ulo Ca l , 



(8) 



where 



U dlpole = -~ / / P(r) • A(r - r') ■ P(r') d 3 rd 3 r' (9) 



and 



Ui oca i — 



^ J |P(r)| 2 d 3 r. 



(10) 



Consider the work done by bringing in from infinity a 
small change in the charge density dp = —divSV. The 
change in Coulomb energy will be given in terms of the 
potential $(r) according to 



SU = 



= J $(r)Sp(r)d 3 r 

= - J $(r)dw<SP(r)d 3 r 

= J grad$(r) • <5P(r)d 3 r. 



(11) 



The electric field E = grad$ may then be computed 
from the functional variation 



SU 



J E(r) ■ <5P(r)d 3 r, 



47T. 



E(r) = F(r) - yP(r). 



(12) 



The local electric field F(r) is due directly to the dipole- 
dipole interaction function A(R) in Eq.Q via the vari- 
ation 



SU, 



dipole 



-J F(r) • <5P(r)d 3 r, 
F(r) = J A(r - r') • P(r')dV. (13) 

The remaining (47r/3)P(r) term on the right hand side of 
Ea. (|12fl is due to the local polarization energy in Ea. (|10fl . 
i.e. 

SUi OC al = SU — SUdipole 

= - f{E(r) -F(r)} • (5P(r)d 3 r 
An f 

= — / P(r) -^P(r)d 3 r. (14) 

Ea.l|Tl)l implies the Clausius-Mossotti local field Eq.JTJ. 

Now suppose that the individual fluid molecules ex- 
hibit a mean electric dipole moment response to a local 
field as described by the polarizability a via 



The dipole moment per unit volume (corresponding to a 
volume v per molecule of fluid) then reads 



„ p aF a / 47T 
P = - = = - E + — P 



V V V 

Thus 

P = XE, 

with a polarization susceptibility 

a 

X ■ 



(16) 



(17) 



(18) 



v i> — (47r/3)a^ 
In terms of the dielectric constant e, 

D = E + 4ttP = £E, (19) 
one finds the Clausius-Mossotti prediction 

'v + (87r/3)a N 



S = 1 + 47TX 



v — (47r/3)o! 



(20) 



All that is needed to calculate e (in this model) is the 
polarizability a for an isolated molecule and the volume 
per molecule v in the fluid. 



III. THERMODYNAMIC STABILITY 

Landau and Lifshitzfl^ prove a non-trivial theorem in 
their treatise on the electrodynamics of continuous me- 
dia. The theorem asserts that the dielectric constant of 
a material obeys the inequality 

1 < e < oo (thermodynamic stability). (21) 

We note (in passing) that the magnetic analogue of the 
theorem (with B = /2H) asserts only the weaker magnetic 
permeability inequality < /i < oo. Thus, diamagnetism 
(0 < /i < 1) can exist but diaelectricity (0 < e < 1) can- 
not exist. If one applies the thermodynamic stability 
criteria of the Landau-Lifshitz Eg. (|2ip to the Clausius- 
Mossotti Ea. H20|) . then it becomes apparent that the re- 
gion of stability for the model may be described by the 
parameter 



V 



{ 47TClA 



r] < 1 (thermodynamic stability). 



(22) 



Debye defines a polar molecule as one whose polariz- 
ability has a temperature dependence (in the dilute gas 
phase) of the form 



3k B T 



(23) 



By plotting a as a function of inverse temperature one 
finds a linear slope yielding the magnitude of a "perma- 
nent" dipole moment 



p = aF. 



(15) 



(24) 
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If the polarizability of a single molecule (say in the gas 
phase) exhibits an appreciable permanent dipole moment 
[/,, then the molecule is called polar. The water molecule 
H2O is an important example of a polar molecule. In the 
dilute gas phase, a water molecule exhibits a temperature 
dependent polarizability of the form in Eq.JSHJ) with 



1.494 x 1(T 24 cm 3 



a {H 2 O) 

H{H 2 0) w 1.855 x 1(T 18 Gauss cm 3 . (25) 

For room temperature and atmospheric pressure, well 
known properties of water [l3j include 



CWer » 28.696 A 3 
V wat er » 30.014 A 3 



(Gas Phase), 
(Liquid Phase), 



(26) 



and 



ywate 



47ra„ 



3v v 



4.0054 > 1 (unstable). (27) 



In TABLE [I] we have listed some liquids along with a 
computation of whether or not the conventional Clausius- 
Mossotti model yields a stable polarization disordered 
state. 

Room temperature (and pressure) water is an impor- 
tant example of a substance which is in an unstable 
regime of the Clausius-Mossotti model. Many other po- 
lar liquids also tend to be in unstable regimes. In their 
studies of polar liquids, the well known phy sical chemists 
OnsagerfiH and (later) KirkwoodflR lit) ] tried to im- 
prove on the Clausius-Mossotti theory in such a manner 
that the instability would be avoided. The attempts were 
not entirely successful. 

More recently, there have been studies which assert 
that polar liquids (such as water) may be found in a 
"superradiant" state wherein ordered domains exist. We 
regard these more recent proposals as likely to be true. 
The superradiant ordered sub-domains of the liquid carry 
a net polarization Pq. From the viewpoint of Clausius- 
Mossotti local fields, the polarized sub-domains should in 
reality be present. In deriving this result below, we shall 



TABLE I: Listed below are some liauids|l3|. The stability 
properties are evaluated on the basis of the Clausius-Mossotti 
model Eq.||55J with 77 = (4ira/3v). 



Chemical 


n 


Clausius-Mossotti 


C2H5OC2H5 


0.5486 


stable 


HCONH2 


6.1498 


unstable 


HCN 


4.7479 


unstable 


HF 


3.6503 


unstable 


HCONHCH3 


5.2093 


unstable 




0.3156 


stable 


PCh 


0.6728 


stable 


H 2 


4.0054 


unstable 


CHCh 


0.5538 


stable 



neglect a® in Eq. l|23|l compared with the thermal term, 
i.e. we employ the Langevin model form 



3k B T 



(Langevin) . 



(28) 



This is certainly a good approximation for water wherein 
aoiH^O) « a water « cut- In what follows we shall also 
review how Ea. (|28l) is derived. 



IV. ORDERED THERMODYNAMIC PHASES 

For a permanent dipole moment p = /in (where n 
is a unit vector), the partition function describing the 
interaction with the local field F is given by 



oP -F/fc B T 



d 2 n 

ill 



(29) 



In Ea. (|29J) . d 2 n is a solid angle about the n direction. 
Explicitly, 



Z = 



^ Fcose/kBT sin6d9, 



yielding a free energy per unit volume of 

'smh(nF/k B T) 



(jiF/k B T) 



(30) 



(31) 



The resulting polarization P = — (0A/dF)T is parallel to 
F and has the magnitude originally derived by Langevin, 
i.e. 



P = ^ 



C ° tH (/S) 



Note that 



Oi T ,. (dP 
= hm — — 

V F^O \ dF 



k B T 



3vk B T 



(32) 



(33) 



in agreement with Ea. l|28l) . 

From Eqs.Q and we may compute (in parametric 
form) the polarization equation of state as follows: 
(i) Define the parameters 



n 
v 



4TTCtT [iF 

-5 — ) x - - — 

6V 
[IP 



k B T 
fxF 



■qk B T k B T 
(ii) Eqs.QJ and now read 



(34) 



z — 'Srj 



coth(z) 



9 

V ~ 47 



coth(z) 

z 



(35) 
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(iii) By eliminating the parameter z in Eas. (|35(l one de- 
fines the function 

y = F(x,ri) (36) 
and thereby the non-linear polarization equation of state 



P 



~9tT \k^f ' 9vk B T 



(37) 



The polarization is plotted in FIG. 1 for room temper- 
ature toluene (which is disordered) and water (which is 
ordered) . 

The polarization equation of state in the Langevin- 
Clausius-Mossotti model is of the typical "mean field" 
variety. There exists a critical temperature 



T r = 



(38) 



P 

c 

»> 0.5 



(T/TJ 

FIG. 2: Shown is the remnant polarization Po(T) for temper- 
atures below the critical temperature. The ordered phase of 
the Langevin-Clausius-Mossotti polar liquid is ferroelectric. 
The ordered net polarized regions should appear in meso- 
scopic domains. The direction of the remnant polarization 
n = (Po/Po) should vary depending upon the domain. 



Above the critical temperature, the susceptibility is pos- 
itive, 

( £ V L )^^(g) T >0fo r T>r c , (39, 

in accordance with the Landau-Lifshitz Ea. i|21[l . Below 
the critical temperature, there is an unphysical "xt < 0" 
portion of an equation of state loop which must be sub- 
ject to an "equal area" construction to deduce the 
proper stable phase. There exists below the critical tem- 
perature a remnant polarization 

P (T) = lim + P(E, T) > for T < T c . (40) 



a. 




(E/E T ) 



FIG. 1: The Clausius-Mossotti "dimensionless" polarization 
(vP/[i) is plotted as a function of the electric field, in units 
of the thermal electric field Et = (feT/(i). Toluene is sta- 
ble with a positive susceptibility xt = (dP / 8E)t ■ Water is 
unstable as indicated by the unphysical "xt < 0" portion of 
the loop in the equation of state. The stable phase of water 
may be constructed using an "equal area" rule. 



The parametric equations for the remnant polarization 
are given by 



T _ 3 
T c z 
vPo = 



coth(z) - 
coth(z) 



(41) 



where the parameter range is < z < oo. The ordered 
phase of the Langevin-Clausius-Mossotti polar liquid is 
ferroelectric with a polarization pointing in a unit vector 
direction n = (P /P ). The magnitude of the remnant 
polarization is plotted in FIG. 2. 

In this picture, water at room temperature contains or- 
dered ferroelectric domains with a polarization pointing 
in the direction n. From one domain to another, the di- 
rection n of the polarization is changing. The situation 
is shown schematically in FIG. 3. Inside of a domain, 
there is a residual local electric field Fo and polarization 
Po related by 



'o = 



47r \ ,> 
Y )Po. 



(42) 



The domain picture has many interesting features only 
one of which we shall discuss here. If one asks what 
will happen to a noble gas impurity (say Argon) when it 
is dissolved in water, then one concludes that the atom 
will be drawn inside the ordered domain. The effective 
potential|l8( on the atom is (roughly) a step function 
attracting the impurity into the ordered domain, 

^Uimpurity (r) — ~^(%impurity |F()(r) | . (43) 

In FIG. 3 we show two of the domains with an impu- 
rity. The notion that an ordered domain is clustered 
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FIG. 3: Shown schematically is a system of polarized domains 
in water. In each domain there exists a local electric field and 
a polarization vector related by Fo = (47r/3)Po. Between 
the ordered phase domains are unpolarized disordered water. 
Two of the domains are shown to contain a noble gas atom 
(such as Argon). The low entropy ordered domain about a 
noble gas impurity has been called an "iceberg" and was un- 
expected when it was first measured thermodynamically. 

around a dissolved noble gas atom (say argon) has con- 
sequences for the thermodynamic impurity equations of 
state. The impurity entropy has in fact been particu- 
larly well analyzed by Frank and Evans p^.l2fj| . Suppose 
that a dilute mixture of Argon in water is such that the 
vapor phase mixture is in thermal equilibrium with the 
dilute liquid mixture. An impurity atom in water va- 
por is in a disordered environment. The impurity atom 
in liquid water will be found in an ordered environment 
within a domain. Impurity entropy is then lost when the 
impurity atom moves from the vapor to the liquid. To 
use the words of Frank and Evans (describing this un- 
expected entropy loss) "the water builds a microscopic 
iceberg around the non-polar molecule" . There is no ev- 
idence that an ordered domain about a noble gas atom 
impurity is a piece of ice. The low entropy polarized do- 
main predicted by the Langevin-Clausius-Mossotti model 
appears to explain the so called "iceberg effect" . 



V. CONCLUSIONS 

The Clausius-Mossotti model together with the el- 
ementary statistical thermodynamics of Debye and 
Langevin yields a ferroelectric phase transition for many 
polar liquids including water. The notion of ferroelec- 
tric domains within the a polar liquid with a net polar- 
ization Po and internal electric field Fo = (4-7r/3)Po is 
clearly implicit in the model and can be derived above in 
a fairly simple and straight forward manner. The phase 
transition may be described using basic concepts from 
elementary physics. 

Nevertheless, reasoning with classical electrodynamics 
and classical statistical thermodynamics yields only par- 
tial insights concerning the complete microscopic theory. 
A detailed theory involves quantum statistical mechan- 
ics and quantum electrodynamics. For example, why 
does a polar H2O molecule have a permanent electric 
dipole moment? If one starts by calculating the quan- 
tum ground state |0) of the molecule employing only the 
non-relativistic Coulomb Hamiltonian, then parity con- 
servation requires the vanishing of the mean dipole mo- 
ment p = (0 |p| 0) = 0. The notion of a permanent 
electric dipole moment requires a break in parity (space 
inversion) symmetry. Such a symmetry break can come 
about due to interactions with the quantum electrody- 
namic field. Similarly, the dipolc-dipole interaction takes 
place (in the quantum electrodynamic theory) due to the 
exchange of virtual photons between two dipoles. 

That the polarized domains have a net internal elec- 
tric field Fo, really means that the exchanged pho- 
tons form a Bose condensate. The Bose condensate is 
also present in a model conventionally used to describe 
superradiance |2l) . In spite of the long time popularity of 
the Clausius-Mossotti model, it was only through a st udy 
of superradiance that the ordered domains were first 
theoretically discovered. It is hoped that the above sim- 
ple Clausius-Mossotti model considerations serve to clar- 
ify the mechanism for ordering in some polar liquids. 
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